A Twisted Kink Crystal in the Chiral Gross-Neveu model 
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We present the detailed properties of a self-consistent crystalline chiral condensate in the massless 
chiral Gross-Neveu model. We show that a suitable ansatz for the Gorkov resolvent reduces the 
functional gap equation, for the inhomogeneous condensate, to a nonlinear Schrodinger equation, 
which is exactly soluble. The general crystalline solution includes as special cases all previously 
known real and complex condensate solutions to the gap equation. Furthermore, the associated 
Bogoliubov-de Gennes equation is also soluble with this inhomogeneous chiral condensate, and the 
exact spectral properties are derived. We find an all-orders expansion of the Ginzburg-Landau 
effective Lagrangian and show how the gap equation is solved order-by-order. 



I. INTRODUCTION 



^ ^ ^ In a recent Letter [1], the authors found a new self-consistent crystalline condensate solution to the gap equation 
^ ' of the massless chiral Gross-Neveu model For this complex chiral condensate, the amplitude is periodic and 
■ the phase winds by a certain angle over each period. Our approach is based on the observation that a carefully 
motivated ansatz for the associated Gorkov resolvent reduces the gap equation to a simple ordinary differential 
equation, an explicitly soluble form of the nonlinear Schrodinger equation (NLSE). In general, the gap equation for 
an inhomogeneous condensate is a highly nontrivial functional differential equation, so the reduction to the NLSE 
represents a significant simplification. This resolvent approach is complementary to the inverse scattering approach 
[3)0|> which also dramatically simplifies the gap equation, but which was not developed for periodic inhomogeneities. 
_ , Our resolvent method is based on an extension, to complex and periodic condensates, of the work of Feinberg and 
Zee Q . The general solution to the nonlinear Schrodinger equation contains all previously known self-consistent 
(— I ' condensates of the massless Gross-Neveu models [both chiral and non-chiral] as special cases: the single real kink [§] , 
I— ' I the single complex kink [1] , the real kink crystal Q , the complex chiral spiral [1| , and also yields a new complex kink 
crystal [l|]. In the language of condensed matter physics, this crystalline condensate is a new solution of the Eilenberger 
^ > equation (for the Gorkov resolvent) @, and we also present here the complete exact solution of the associated 
Qv^ Bogoliubov-de Gennes [l^ equation H%l) = for this system. The Eilenberger and Bogoliubov-de Gennes equations 
are fundamental elements of the treatment of a wide class of interacting fermion systems, which are important in 
many branches of physics, ranging from particle physics, to solid state and atomic physics [HI, [H, [l^ [3| ■ Important 
C^l [ paradigms include the Peierls-Frohlich model of conduction ,15], the Gorkov-Bogoliubov-de Gennes approach to 
\^ ■ superconductivity [lo|, and the Nambu-Jona Lasinio (NJL) model of symmetry breaking in particle physics [l6l|. Here 
we study a 1 + 1-dimensional version of the NJL model, the NJL2 model [also known as the chiral Gross-Neveu model, 
00 , XGN2]. This model has been widely studied as it exhibits asymptotic freedom, dynamical mass generation, and chiral 

' symmetry breaking 0, H, 0, @ . 
^ , Our primary physical motivation for studying the gap equation of the massless chiral Gross-Neveu model, the 
NJL2 model, is to understand the (T, ^) phase diagram of this system. Somewhat surprisingly, the phase diagram 
of this system is not yet fully understood. A gap equation analysis based on a homogeneous condensate suggests 
its phase diagram is the same as its discrete-chiral cousin, the original Gross-Neveu (GN2) model 0], while more 
recent work finds an inhomogeneous Larkin-Ovchinikov-Fuldc-Ferrell (LOFF helical complex condensate ( "chiral 
spiral" ) below a critical temperature [11 . In [l| , a Ginzburg-Landau approach was used to show that in a region of 
the phase diagram the free energy is lower for a complex kink crystal, compared to a uniform condensate or a chiral 
spiral condensate. In this paper we present the details of the complex crystalline condensate of the NJL2 system, and 
also the exact spectral properties of fermions in the presence of such a crystalline condensate. This information will 
subsequently be used to study the free energy exactly, without resorting to the Ginzburg-Landau approximation. 

This state of affairs should be compared and contrasted with the case of the original Gross-Neveu model 0, to 
which we refer as the GN2 model, which has a discrete chiral symmetry rather than the continuous chiral symmetry 
of the NJL2 model. In the GN2 model, the phase diagram has only relatively recently been solved in the particle 
physics literature, analytically and exactly by a Hartree-Fock analysis and numerically on the lattice [3]. There is a 
crystalline phase at low temperature and high density, and this phase is characterized by a periodically inhomogeneous 
(real) condensate that solves exactly the gap equation. This phase is not seen in the old phase diagram which was 
based on a uniform condensate (lol . [20| . Interestingly, some hints of a problem with the homogeneous condensate 
assumption were found already in an early lattice study [21]. This discrete-chiral GN2 model (with vanishing bare 
fermion mass) turns out to be mathematically equivalent to several models in condensed matter physics 0]: the real 
periodic condensate may be identified with a polaron crystal in conducting polymers [ill, US with a periodic 
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pair potential in quasi ID superconductors [2J, |25|, |2g| , and with the real order parameter for superconductors in 
a ferromagnetic field [2 7| . The Gross-Neveu models also serve as paradigms of the phenomenon of fcrmion number 
fractionalization [H, l29l[30l. [3lj | . 

Here, we consider the massless chiral Gross-Neveu, or NJL2, model in 1 + 1 dimensions with Lagrangian 0, 



(1.1) 



This system has a continuous chiral symmetry under ip ~~f e*^ "ip. We have suppressed summation over N flavors, 
which makes the semiclassical gap equation analysis exact in the N 00 limit, a limit in which we can consistently 
discuss chiral symmetry breaking in 2D [s^. l33l|. The original Gross-Neveu model, the GN2 model without the 
pseudoscalar interaction term (^'tpi'j^tp'^ , has a discrete chiral symmetry ip — > j^ip. 

There are two equivalent ways to find self-consistent static condensates. First, introduce bosonic scalar and pseu- 
doscalar condensate fields, a and tt, which we combine into a complex condensate field, defined either through its real 
and imaginary parts, or via its amplitude and phase: 



Integrating out the fcrmion fields we obtain an effective action for the condensate A as 



5eff = / l^l'-*^lndet 

The corresponding (complex) gap equation is 



5A{x)* 



In det 



^(l-7^)A-i(l+7^)A' 



(l-7^)A(:r)--(l-t-7^) A* {x) 



(1.2) 



(1.3) 



(1.4) 



If the condensate is constant, as is usually assumed, it is straightforward to evaluate the determinant and solve the 
gap equation [l^, [l^, [23] ■ When the condensate is inhomogeneous this is a much more difficult problem. Dashen, 
Hasslacher and Neveu [3] used inverse scattering to find kink-like static but spatially inhomogeneous condensates for 
the gap equation of the GN2 model (where there is no pseudoscalar condensate, so A is real) . Shei Q extended this 
inverse scattering analysis to the chiral Gross-Neveu model, the NJL2 model, and found a spatially inhomogeneous 
complex kink. A new approach to the inhomogeneous gap equation, based on the resolvent, was developed by Feinberg 
and Zee 0] and applied to the kink solutions of both the GN2 and NJL2 models. For the GN2 model, Thies used a 
Hartree-Fock approach to find a periodic extension of the real kink solution, motivated by analogous inhomogeneous 
condensates in condensed matter systems [7| . In py] , the present authors showed that the complex gap equation (jl.4p 
can be reduced in an elementary manner to a soluble form of the nonlinear Schrodinger equation. The general solution 
contains all previously known inhomogeneous condensates (real and complex), and yields a new crystalline extension 
of Shei's complex kink. 

A second approach to finding a self-consistent condensate is to solve the relativistic Hartree-Fock problem Hijj — Eip, 
with single-particle Hamiltonian 



i/ = -i7 — +7 
ax 



l(l-7^)A(x) + i(l + 7') A*(x) 



and subject to the consistency condition 



(1.5) 



(1.6) 



We choose Dirac matrices 7*^ = ci, 7^ = —i<J2, 7^ — f^s, to emphasize the natural complex condensate combination 
in p.2p . Then the single-particle Hamiltonian is 



H 



A(x) 



"' dx 



(1.7) 



This Hamiltonian is also known as the Bogoliubov-de Gennes (BdG) Hamiltonian, and we will refer to the associated 
spectral equation 



Hip = 



(1.8) 
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as the Bogoliubov-de Gennes (BdG) equation. 

In Section |TT] we review the reduction of the functional gap equation to the nonhnear Schrodinger equation (NLSE), 
and in Sections IIIII and ITVl we present the real and complex condensates obtained from solving the NLSE. In Section 
IVl we show that the associated Bogoliubov-de Gennes equation can also be solved exactly, and we derive the exact 
single particle spectrum and density of states. In Section IVll we verify the consistency of our solutions by solving the 
gap equation in the Hartree-Fock approach. An all-orders Ginzburg-Landau expansion of the free energy is presented 
in Section IVIIl and we show that the inhomogeneous gap equation is satisfied order-by-order in an interesting and 
nontrivial way. In a concluding section we review our results and discuss implications for the phase diagram of the 
chiral Gross-Neveu model. 



II. REDUCTION OF FUNCTIONAL GAP EQUATION TO NONLINEAR SCHRODINGER EQUATION 

In this Section we review the reduction '1] of the functional gap equation (jl.4[) to the nonlinear Schrodinger equation. 
The key quantity in our approach is the coincident limit of Gor'kov Green's function, or the the "diagonal resolvent" : 

Rix;E)^{x\j^\x) . (2.1) 

The resolvent (I2.1|) is clearly a 2x 2 matrix. For a static (but possibly spatially inhomogeneous) condensate, all spectral 
information is encoded in the resolvent. Indeed, the spectral function characterizing the single-particle spectrum of 
fermions in the presence of the condensate A (a;) is 

p(E) = -Im Ttd,x [R{x; E + ie)] , (2.2) 

TT 

where the trace is a Dirac trace as well as a spatial trace. 

Our first, very simple, observation is that the form of the BdG equation (jl.Sp places very strong constraints on 
the possible form of R{x\ E). For any static condensate A(a;), R{x\ E) must satisfy the following algebraic conditions 
[these are explained in more detail in Appendix A]: 

R ^ r!^ (2.3) 
tr D{R{x;E)a^) - (2.4) 

deiR{x;E) = . (2.5) 

Furthermore, R{x] E) must satisfy the first-order differential equation 



d 

—R{x;E)(T3 = I 



E -A(a;) 
A*(a;) -E 



, R{x; E) CT3 



(2.6) 



In superconductivity, (|2.6p is known as the Eilenberger equation [sH] , and in mathematical physics as the Dik'ii 
equation l341|. These conditions (|2.3p - (|2.5p . and the Eilenberger equation l|2.6p . all follow from the simple fact 
[3^ . [35I l36l| that for the one-dimensional BdG equation, which involves derivatives with respect to the single variable 
X, the Green's function can be expressed as a product of two independent solutions to ([TT 



R{^; ^) = ^ {^^^2 + ^2^-0 (2.7) 

where W is the Wronskian of two independent solutions 7/^1.2: W — i4'i<^2'4'2- 

The next step is to note that the gap equation provides further information about the possible form of the resolvent, 
and this is enough to motivate a specific ansatz form There are two ways of viewing the gap equation (|1.4p in 
terms of the resolvent. First, for a static condensate we can write the log det term in the effective action (|1.3p as 
minus the grand canonical potential, in terms of the single-particle spectral hmction p{E): 

r dEp{E)ln(l + e-^'^^-^'A . (2.8) 
P J —00 



All dependence on A(a;) resides in the spectral function p{E), via (|2.2p . Therefore, inserting this into the gap equation 
p.4p . this relates A{x) to the diagonal entries of R{x;E). Further, as a consequence of the condition (|2.4p . these 
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diagonal entries are equal. So, the simplest natural solution to the gap equation is for the diagonal entries of R{x; E) 
to be linear in |A(x)p. A second way to view the gap equation is to evaluate the functional derivative in (jl.4p . which 
for a static condensate leads to: 

A(a;) = -iNghro^E [7° (l + 7'') R{x; E)] (2.9) 

The Dirac trace then relates the off-diagonal entries of R{x;E) to A(x). Since R is hermitean, these off-diagonal 
entries are complex conjugates of one another. 

Summarizing, R{x; E) must be a hermitean 2x2 matrix with equal diagonal entries, such that [after the spatial 
and energy trace] the variation of the diagonal terms is proportional to A(a::), and with off-diagonal terms linear in 
A(x), after the energy trace. This suggests taking the resolvent to be of the form 

R{x,E) -f^[E)y^ 6(£;)A*(:e) a{E) + \^{x)\'' ) ^^'^"^ 

where a{E), b{E) and N{E) are functions of i?, independent of x, and are to be determined. However, this ansatz 
cannot describe inhomogeneous condensates because the only solution of this form consistent with (|2.5p is a condensate 
with constant magnitude, independent of x. Indeed, taking A to be constant (and by a global chiral rotation, real), 
A = M, the solution to is simply 

^A i i^l (2-11) 



2VAf2 - £2 \M E 

as is familiar. This example also illustrates that the hermiticity condition (|2.3p must of course be interpreted with 
the appropriate ie prescription for the energy. 

To find inhomogeneous condensates, we suggested in [l| to extend the ansatz (|2.10p to include a first derivative 
term in the off-diagonal entries: 

This is the simplest extension of ()2.10p that is consistent with the various algebraic constraints and with the Eilenberger 
equation (|2.6p . Indeed, substituting the ansatz (|2.12p into the Eilenberger equation (|2.6p . we see that the diagonal 
entry of this equation is identically satisfied, while the off-diagonal entry implies that A(a;) must satisfy the following 
nonlinear Schrodinger equation (NLSE) [and its complex conjugate]: 

A"-2|A|2A-|-i(6-2£;)A'-2(a-£;6)A = . (2.13) 

Two comments are in order. First, it is not immediately obvious that R{x] E) in (|2.12p can satisfy the normalization 
condition (|2.5p for an inhomogeneous condensate, since 

detR{x;E)=^{\A\^-\A'\^ + {2a-b^)\A\^+ib(A'A*-AA*')+a''} (2.14) 

Remarkably, the NLSE (|2.13p implies that (det i?(x; £^)) is constant: 

d fdetR{x;E)\ ^ (2| A|2 + 2a - 6^) (| A|2)' - (A"A*' A'A*") + ife (A"A* - A*"A) 

= (2.15) 



dx V AA2 



where we have used the fact that the NLSE ([2^3]) implies that (A"A*' + A'A*") = (2|A|2 + 2a - 2Eb) (|A|2)', and 

that (A"A* — A*" A) — —i{b — 2E) (|A|2)'. Since deti?(a::;£') is constant, the normalization in (12. 5p can be achieved 
by suitable choice of N{E). Second, while the ansatz (|2.12p automatically satisfies the x dependence of the gap 
equation in its form coming from (|2.8p [because the trace of R is, by construction, linear in |A|2], it doesn't satisfy 
the other form of the gap equation (|2.9p . until the energy trace is performed. This is because of the A'(a;) terms in 
the off-diagonal. In Section IVll we show that this form of the gap equation is indeed satisfied because the coefficient 
of the A'{x) term vanishes due to the energy trace. 

Thus, we have reduced the very difficult problem of solving the functional gap equation (|1.4p for a self-consistent 
condensate A{x) to the much simpler problem of solving the NLSE for A(x). In fact, the NLSE (I2.13P is explicitly 
soluble, as is discussed in the following sections, in which we describe first the real solutions [relevant for the GN2 
model] , and then the complex solutions [relevant for the N JL2 model] . 
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III. REAL SOLUTIONS OF THE NLSE 

In this Section we recah the previously known real solutions to the gap equation, and show how they fit in with the 
NLSE (Pr^ and the resolvent form in (ITT^ . 

A. Homogeneous condensate 

If the condensate is constant, then by a global chiral rotation it can be taken to be real: 

A{x) = m (3.1) 

This clearly satisfies the NLSE ([^1^ . and we find 

a{E) = 2E^ - 
h{E) = 2E 

The spectrum of the associated BdG equation p.8p is that of a free fermion with mass m, with positive and negative 
energy continua starting at E = ±m, the mass scale being set by the amplitude of the condensate. 

B. Single real kink condensate 

A well known nontrivial solution to the gap equation is the single (real) kink 0]: 

A.{x) = m tanh(mx) (3-3) 

This satisfies the NLSE 

A" - 2 + 2m^A , (3.4) 
and so we deduce the exact diagonal resolvent to be of the form p.l2|) with 



'2 



a{E) = 2E 
h{E) = 2E 

The spectrum of the associated BdG equation (jl.Sp has positive and negative energy continua starting aX E — ±m, 
together with a single bound state located at i? = 0, at the center of the gap. This mid-gap zero mode has many 
important consequences in a variety of branches of physics (28l . [30| . 

C. Real kink crystal condensate 

A periodic array of these real kinks also provides a solution to the gap equation. This solution describes a polaron 
crystal in polymer physics (2^ . [23| , a periodic pair potential in inhomogeneous superconductors 0, [2^, [l^l , and the 
crystalline phase of the Gross-Neveu model . Define (the peculiar looking scaling will become clear below) 

r- 2r7i / 2m \ 
A(x)=V^^^sn(^^^.x;.j (3.6) 

where sn is the Jacobi elliptic function [33, [H, [3^, with real elliptic parameter < < 1 . The sn function has 
period 2Kiv), where K(i^) = J^^'^ (l — ly sin^ t) dt. When v ~ \ (|3.6p reduces to the single kink condensate in 



3.3p . The periodic condensate (|3.6p satisfies the NLSE 



2A3 + (l + i.) ( ^I^Va-O . (3.7) 



FIG. 1: The real kink crystal condensate (|3.6|l plotted for elliptic parameter v — 0.99 [solid, red curve], and for v = 0.1 [dashed, 
blue curve]. For small the condensate has the Larkin-Ovchinikov-Fulde-Ferrell (LOFF) form of a small amplitude sinusoidal 
condensate, while for — > 1 the condensate resembles an array of kinks and anti-kinks. 



Thus, we deduce the exact diagonal resolvent to be of the form (I2.12p with 

a{E) - 2E'-2m \ (3.8) 

(1 + Vi^)2 

b{E) = 2E (3.9) 

^E) - \ J _ (3.10) 

4 Vm2 - E^ ^E^ - m2(i^)2 

This periodic condensate is plotted in Figure [TJ Note that over the period x G |-_ k(i/)(i+v^) ^ k.(v){i+^) -^^ 
condensate is shaped like a single kink. This reflects the expansion of the Jacobi sn function in terms of an array of 
periodically displaced tanh functions: 

oo 

sn(x;z.) = — 1^ J2 {-l)-tanh(^^{x-2nK)) (3.11) 

^ n— — oo 



where we use the standard notation K'(z/) = K(l — z^). In the infinite period limit {v —^ 1), the interval 

j_ , ^^'"''2m~^'' ] i^s-PS to the whole real line, and K' 7r/2, and so the kink crystal (|3.6p reduces pre- 

cisely to the single kink condensate in p.3p . 

It is worth noting that this periodic kink crystal (|3.6p can be written in an equivalent, but different looking, form, 
by use of a Landcn transformation 37, 38, 39, 40] of the Jacobi functions. That is, by rescaling the elliptic parameter 
v together with the argument mx, we can write 

^ sn{mx;i>)cn{mx;D) , iy/i^ 
Aixj^mv — — ; ^ j-, i=rT^ ■ (3.12) 

^ ^ dn{mx;v) (1 + ^)2 ^ ' 

This is the form in which this periodic kink solution is presented in the work of Thies et al on the crystalline phase 
of the Gross-Neveu model, while the form (|3.6p was used in the condensed matter literature in [l^, [H, 0, [2^, [la, . 
The spectrum of the associated BdG equation (|1.8p has positive and negative energy continua starting at i? = ±m, 

together with a single bound band in the middle of the gap, with band edges at i? = ± (^yT^) "^^^^ band lies 
symmetrically in the center of the gap. The spectrum is plotted in Figure [2] as a function of the elliptic parameter v. 
Notice that there is just one bound band in the energy gap, and when v ^ 1 [the infinite period limit], the bound 
band at the center of the gap contracts smoothly to the single bound zero mode of the kink condensate. 



FIG. 2; The band spectrum of the real kink crystal, showing the positive and negative energy continua and the bound band, 
as a function of the elliptic parameter v. The energy is given in units of the scale m. The infinite period limit is — > 1, where 
the bound band shrinks to a single bound level at i5 = 0, the familiar zero mode of the kink condensate. 



IV. COMPLEX SOLUTIONS OF THE NLSE 
A. Single plane wave condensate 

The simplest complex solution to the NLSE is a single plane wave: 

A = me*'^ (4.1) 

This satisfies the NLSE ((2T3)l with 

a{E) = 2(i?-|)'-m2 
h{E) = 2E-2q 

N{E) = \ , ^ = (4.2) 

{E - q/2)^m'^ ~ {E ~ q/2f 

This plane wave condensate behaves just like a constant one, but with the energy shifted by g/2, as can be seen by 
making a local chiral rotation: 

A^e-*«^A , ^^e-'^V^Ts^ _ (43) 

It is clear from the BdG equation (|1.8[) that such a transformation has the effect of shifting the entire energy spectrum 
by q/2. This illustrates an important point: for complex solutions A(a::) of the NLSE, one can always multiply by an 
arbitrary plane wave phase factor e*''^, and this simply corresponds to shifting the entire energy spectrum. In general, 
a local chiral rotation through angle a{x) leads to a local chemical potential [4l|, [i^l, or local electrostatic potential, 
Aq{x) = ia'(a;), and therefore a local electric field E{x) = —^a"{x). For the single plane wave condensate in (|4.ip . 
a{x) is linear in x, and so there is no associated electric field. 



B. Single complex kink condensate 



Shei [4| found a solution to the gap equation for the NJL2 model, in which both the scalar and pseudoscalar 
condensates have a kink-like form: 

a{x) = m [cos^(6'/2) + sin^(e'/2) tanh (m sin (9/2) x)] 

tt{x) = -y sin(6') [1 - tanh(TO sin(6'/2) a;)] (4.4) 
where 9 G [0, 27r] is a parameter. These kinks are plotted in Figure[3l This complex kink (|4.4|) has also been extensively 
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FIG. 3: Plots of the real and imaginary parts (j(x) and tt{x) of the complex kink condensate in (|4.4p for three different values 
of the winding parameter 0. The scalar condensate cr{x) [solid, red curves] winds from m at x = — oo, to mcos(6') at a; = +oo, 
while the pseudoscalar kink n{x) [dashed, blue, curves] winds from —msm{6) to as a: ranges from x = —oo to x — +oo. The 
plots are for 9 — it/10, 9 = 7r/2 and 9 — Qtt/IO, and a and vr are plotted in units of m. 



studied in the resolvent approach by Feinberg and Zee 5]. In our analysis it is more natural to combine these into 
the complex condensate A{x) — a{x) — Z7r(a;), [as in ()1.2|1 ]: 



, cosh (m sin (6'/2) x — i6'/2) .n/n 
^ ' cosh(m sin(6'/2) x) 



(4.5) 



This complex form is plotted in Figure [H This illustrates the role of the parameter G G [0, 27r] as the net rotation 
angle of the kink as x goes from — oo to +oo: 

A(a; = +CX)) = e"*''A(a; = -oo) (4.6) 

Observe that when 6 — n, the complex kink (|4.5p is in fact real, and reduces to the familiar real kink solution in 
p.3|) : this real kink changes its sign [i.e., rotates through tt] in passing from x — —oo to x — +oo. Another useful 
representation of this kink is in terms of the magnitude and phase A{x) — M(a;)e*'^*^^^: 

M^{x) = [l - sin^ {6/2) sech^ (to sin (9/2) x)] 



x(x) — arctan 



smt 



cos 6 + e2mxsin{e/2) 

The complex kink condensate (|4.5p satisfies the NLSE: 

A" - 2 I A|2 A - 2 i m cos (9/2) A' + 2m^A = 

From this NLSE, we deduce the exact diagonal resolvent to be of the form (|2.12p with 

a{E) = 2E^-2mcos{9/2)E-m^ 
b{E) = 2E~mcos{9/2) 

Af{E) = ^ 



(4.7) 
(4.8) 

(4.9) 



4:{E - TO COS {9/2))Vm^ - E^ 



(4.10) 



FIG. 4; Plot of the complex kink condensate (|4.5p . for 9 = 3n/2, illustrating how the kink winds around zero without the 
amplitude vanishing. The kink is the solid [red] line, and the surface is shown simply to illustrate that both the amplitude and 
the phase are changing. 



M, X 
■■%.5 - 
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FIG. 5: Plots of the amplitude M and phase x of the complex kink condensate (|4.5p . for three different values of the winding 
parameter 9. The condensate amplitude M{x) [solid, red curves] approaches m at x = ±oo, and equals mcos{9/2) at the kink 
center x = 0. The phase x{^) [dashed, blue, curves] winds from S to as a: ranges from x = —oo to x = +oo. The plots are 
for 9 = tt/IO, 9 = 7v/2 and 8 = 97r/10. 



The spectrum of the associated BdG equation p.Sp has positive and negative energy continua starting a.t E — ±m, 
together with a single bound state located a.t E = m cos (0/2). When 9 — tt, where this complex kink reduces to the 
standard real kink, the bound state is once again a zero mode. But for other values of 9 the single bound state lies 
asymmetrically inside the gap, as plotted in Figure [6]. As 9 goes from to 27r, one state moves from the positive to 
the negative energy continuum. As is clear from the previous subsection, we can always multiply the complex kink 
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FIG. 6: Plot of the fermion single-particle spectrum for the single complex kink (|4.5|l . as a function of the winding parameter 
9. Note that for 6 = ir [when the condensate is real] the bound state is at 0, but for all other values of 6 the bound state lies 
asymmetrically in the gap. 



solution (|4.5p by a plane- wave factor e"'^, which has the net effect of displacing the fermion spectrum by q/2, with 
the corresponding simple modifications to the resolvent functions a(E), b{E) and M{E) in (|4.10p . 

At this stage we have shown that Shei's complex kink condensate (|4.5p solves the NLSE, and we have found the 
corresponding exact diagonal resolvent (|2.12p with a{E), b{E) and A/'^E) given in (|4.10p . This agrees with the spectral 
properties derived from inverse scattering [J]- Shei further showed [J| that this complex condensate solves the gap 
equation provided a further restriction is applied to the winding parameter 9. This condition states that 6/{2-k) is 
equal to the filling fraction in the large flavor limit, of the single bound state in the gap by n flavors, with flxed 
as iV ^ oo 



27r 



n 
N 



(4.11) 



In Section IVII we show that in our approach this same condition arises from demanding that the coefficient of the 
A'(a;) term in (|2.9p vanishes after the energy trace, a necessary requirement to satisfy the gap equation. 



C. Complex kink crystal condensate 

A new complex condensate was presented in [ij. This new solution is a periodic array of Shei's complex kink (j4.5p . 
Physically, it is associated with a crystalline phase of the NJL2 system [l], just as the real kink crystal condensate 
(13. 6p is associated with a crystalline phase of the GN2 system 7]. Up to a plane wave factor [as in Section FlV Aj . 
this complex crystalline condensate is the general solution to the NLSE (|2.13p . and all other solutions [both real and 
complex] can be obtained from it by suitable choices of parameters. This solution can be written [6J] in terms of 
Weierstrass elliptic functions [or, alternatively but equivalently, in terms of Jacobi theta functions] 



A(x) = -A 



a{Ax + iK' -i9/2) 



exp [iAx{~i C{i9/2) + i ns{i0/2)) + i 6m/2] 



'a{Ax + iK')a{ie/2) 
The parameter A sets the scale of the condensate and its length scale: 

A = A{m, e, v) EE ~2i m sc {iO/A) nd [iO/A) 



(4.12) 



(4.13) 



where sc=sn/cn and nd=l/dn are Jacobi elliptic functions [37j, |38|. The functions a and C, are the Weierstrass sigma 
and zeta functions [13, [H, [H, H^] , some relevant properties of which are given in Appendix B. We have chosen real 
and imaginary half-periods: wi = K(i/), and W3 = iK' = iK(l — v). Both periods are therefore controlled by the 
single [real] elliptic parameter < < 1. Also, 7/3 = C(*K') is purely imaginary. The parameter 9 € [0,4K'(i^)] is 
related to the angle through which the condensate rotates in one period L = 2K/A : 



A(a; + i) ^ e2*'^A(x) 



where the angle is a function of 9 and v 



(p = K(-*C(*0/2) + zns(*0/2)-^ 



(4.14) 



(4.15) 
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FIG. 7: Plot of the complex kink crystal condensate (|4.12|) . for v = 0.8 and d — 3K(0.2)/2, illustrating how the kink winds 
around zero each period, without the amplitude vanishing. The kink is the solid [red] line, and the surface is shown simply to 
illustrate that both the amplitude and the phase are changing over each period. 




FIG. 8: Plots of the amplitude M{x) and phase x{^) of ths complex kink crystal condensate (|4.12|l . for 9 = 1.6 and f = 0.95. 



Here we used the quasi-periodicity property (|11.2[) of the a function. Note that (p and ij = C(K) are real, and when 
v —^ I, we have (p — *■ —0/2. This crystaUine complex kink is plotted in Figure [7] showing the winding of the kink over 
a period. 

It is also useful to visualize the condensate (|4.12p in terms of its amplitude and phase: A = Me*-^. The modulus 
squared is a bounded periodic function, with period 2K./A: 



|A(x)|2 = {V {Ax + iK') - V (ie/2)) 



(4.16) 



Here we used the quasi-periodicity property pi.2p of the a function, together with the product identity (jll.l2|) 
relating the a and V functions. The phase x can be expressed as 



/.,^x f<j(Ax + iK' + id/2)\ n^e 

= AH C W2) + . ns(../2)). + - m j + \- 



(4.17) 



The amplitude and phase are plotted in Figure [SI Note that the amplitude is periodic while the phase changes by 2if 
over each period. 

The complex crystalline condensate in ()4.12p satisfies the NLSE: 



A" - 2 |A|2 A - i {2Ains{ie/2)) A' - A^ {iV{ie/2) - xis^{ie/2)) A = 



(4.18) 



Comparing this equation with the NLSE (|2.13[) we can extract the functions a{E), b{E) and Af{E) appearing in 
p.l2p . thereby determining the exact diagonal resolvent. To express these functions in a compact form, we define 
some properties of the associated fermionic spectrum for the BdG equation p.Sp . This spectrum has positive and 
negative energy continua starting at i? = ±m, together with a single bound band in the gap, as depicted in Figure [HI 
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FIG. 9: Plots of the single-particle fermion spectrum for the complex kink crystal condensate (|4.12p . for = 0.1 [first plot] and 
v = 0.9 [second plot], as a function of the winding parameter 9. Note that for 9 = 2K' [when the condensate is real] the band 
is centered symmetrically about E = G, but for all other values of 9 the band lies asymmetrically in the gap. In the infinite 
period limit, v ^ 1, the bound band shrinks to a single bound state, and its 9 dependence reduces to that depicted in Figure 
|6]for the single complex kink condensate. 



In contrast to the case for the real kink crystal in Section IIIIl here the bound band is not centered in the middle of 
the gap, but is displaced from the center. The parameter characterizes this asymmetry in the spectrum. The band 
edges are functions of both the winding angle 9 and the elliptic parameter v: 

El = —TO 

E2 = m(-l + 2nc2(i6l/4;i/)) 

^3 = m{-l + 2TLd^{ie/A-v)) 

E4 = +TO (4.19) 

In the infinite period limit [v 1], the band contracts to a single bound state, with E2 = E3 = tocos {0/2), and this is 
precisely the bound state of the single complex kink, as shown in Figure[6] At finite period, but when 6 — 2K'(i/), we 

find E2 = —Ej, = — ^ 1^^ ^ TTi, and the band is centered symmetrically about 0; this is precisely the band spectrum 
of the real kink array in Section IIIIl Thus, we can roughly think of the parameter as setting the offset location 
of the band inside the gap, while the parameter u, together with 0, plays the role of determining the period of the 
crystal, and hence the width of the band. In terms of the band edges, the resolvent functions a(£'), h{E) and M{E) 
that appear in the resolvent (|2.12p take the following simple form 

a{E) = 2E^ - {E2 + E^)E - ~ - 

b{E) = 2E-(E2+E3) 

N{E) = — ^ = (4.20) 

A,/{m'^-E^){E-E2)(E~Es) 



Here the normalization Af{E) is fixed by the property det i? = — 7 



4' 



1. Solution of NLSE by separation into amplitude and phase 

We now comment on two ways to derive this solution to the NLSE (|2.13|) . The first is to separate the NLSE into two 
equations, one for the amplitude and one for the phase. Writing A(.t) = M{x)e'^'-''\ and considering (A*A"- A"*A), 
we immediately find that the phase is related to the amplitude by: 

X»^-l(i,-2i=) + -^ (4,21) 

where Ci is a constant. Note that this is indeed true for all the complex condensate solutions discussed above. For the 
crystalline complex solution in (|4.16l I4.17P the relation (|4.2ip can be verified using the Weierstrass function properties 
listed in Appendix B, in particular equation pi.lSp . 



13 



Next, considering (A* A" + A"* A), we find the following nonlinear equation for M"^: 

M^yy = AM^ + 2 ^4(a -Eb)-^{b- 2Ef^ + C-zM^ ~ ACf (4.22) 

where C2 is another constant. From the form of this equation, comparing with the equation (|11.10p for the Weierstrass 
V function (stI . IssI . IsqI . l40l | , we recognize the solution to be of the form 

= A^[V{Ax + i'K')-V{C:i)] (4.23) 

The shift by iK.' ensures that M'^{x) is bounded for x on the real axis. This explains the result in (|4.16[) . identifying 
C3 = iO/2. The constants Ci and C2 are related to C3 as 

Ci = -'-A^V'{iei2) 

C2 = 2A'^V"{i9/2) (4.24) 
Given this solution for the amplitude, the phase follows from ()4.21|) . using the integration formula [stI ]: 

' '"<^=> ^^,.l„(£<4£±i^V2(^« + '-KK(C3) (425) 



A J ViAx + iK)-ViC3) \(T{Ax + iK + C3) 

Notice that the phase can always be shifted by qx, which amounts to multiplying the solution by a plane wave factor, 
as discussed in Section IIV Al In our solution (|4.12p , the plane wave factor has been chosen so that the associated 
fermion spectrum has the form shown in Figure [9l An additional plane wave factor in A(x) displaces this entire 
spectrum by a constant. 



2. Solution of NLSE as periodic array of complex kinks 

Another way to derive the general periodic complex solution to the NLSE is to make an educated guess for a periodic 
array of Shei's complex kink solution (j4.5p . using known properties of the Weierstrass functions. Shei's complex kink 
condensate (j4.5|) can be written as 

, sinh(TOsin(6'/2)a; + i7r/2-i6'/2) 
= - sinh(^sin(Ux + W2) ^'-''^ 

This can be made quasi-periodic along the real x axis by generalizing the hyperbolic sine function to its doubly-periodic 
form, which is the Weierstrass sigma function. Thus, we are led to try the form 

A(.) = ^e--/^^i^^L±^;:i^e- (4.27) 
a (Ax + iK') a (tO/2) 

■\ . . . 16*773/2 

where we have included a possible plane wave factor , to be determined, and the normalization factors ^^^^^^^ 

have been chosen for convenience. 

Given this form of A (a;), since C(a^) = d/dxlTia{x), we find 

A' = A {({Ax + iK' ~ i6/2) - C,{Ax + iK!) + iX/A] A (4.28) 

and since V = — C'l 

A" = A^ \^[-ViAx + iK' - id/2) + V{Ax + iK')] + [C {Ax + iK' - i6/2) -C{Ax + iK') + i\/Af^ A (4.29) 

Furthermore, using the quasiperiodicity properties of the Weierstrass functions (|11.2|) and the Weierstrass function 
product formula pi.l2p . it follows that 

2 _ 2 {Ax + iK' - ^61/2) cr {Ax + iK' + ie/2) 
' ^"^'^ ^ a2 {Ax + iK')a{i6/2)a{-id/2) 

= A^[V{Ax^iK')~-V{ie/2)] (4.30) 
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Therefore, 



A"-2|ApA ^ yl2{-[p(Aa; + iK'-i6l/2)+P(Ax + iK')] 



+ [C [,Ax + iK' . 
+2P(i6l/2)}A 



i6i/2)-C(^a; + iK') + iA/A]' 



(4.31) 



where note that the relative sign between the first V functions has been flipped by the subtraction of 2|ApA. 
Now we use the remarkable identity (|11.14p that relates the V function to squares of the C function, to find that 



A" - 2A2A = 



iV{ie/2) + {\/A + iC[ie/2)y 



2iA [A/.4 + <(?;6i/2)] A' 



(4.32) 



Thus, A(x) of the form in (|4.27|) does indeed satisfy the NLSE ()2.13p . and one simply needs to match the constants 
algebraically in order to express a{E), b{E) and Af{E) in terms of 9. This determines A in (|4.27p up to an additive 
constant, which can be fixed by matching the fermion spectrum to the form in (|4.19p . 



D. Reduction of general solution to special cases 

To conclude this section describing our new complex crystalline solution ()4.12p . we note that it incorporates all 
previously known solutions as special cases. 



1. Reduction to real kink crystal condensate 

When 9 = 2K', the scale factor A 2m/ (1 + y/T'), and the condensate reduces to the real kink crystal in (|3.6 

2m\fv ( 2m x 
A -> ^sn 



^ sn (m x; i>) en (m X] i>) _ 4^^ 



m V ■ 



ATi{mx\v) ' (1 + v^)2 



(4.33) 



In this limit, the amplitude vanishes at x = 0, and one sees that the kink "winds" through angle 2iy9 — vr by passing 
through zero. For other values of 0, the complex kink crystal winds around zero, but without the amplitude actually 
vanishing. For = 2K', the bound band becomes symmetric about E — 0, because the band edges (|4.19p reduce to: 

Accordingly, the resolvent functions a{E), K^) JV{E) in (|4.20p reduce smoothly when 6 = 2K' to the corre- 
sponding expressions for the real kink crystal in p.lOp . 



2. Reduction to single complex kink condensate 

For general 6, when v ^ I, the complex crystalline condensate in (|4.12p reduces to Shei's complex kink solution 
()4.5p . To see this, first observe that in this limit, A{m,9,v) m sm{6/2), and K diverges, while K' — > tt/2. Thus 
the period 2'K/A diverges. Also, the band edges contract as E2 ^ E^ = mcos{9/2), so the band shrinks to a single 
bound state, whose 9 dependence matches that of the bound state for Shei's complex kink condensate. Furthermore, 
as — !■ 1, the Weierstrass functions simplify: 

a{x) ~* smh{x)e^^^^^ 
({x) — ^ + coth(a;) 



'Pix) -> - + . , 2. , (4.35) 



1 _^ 1 
3 sinh^ (x) 

These relations show that on the interval [— K/A, K/A], in the limit — > 1, the complex crystalline condensate (|4.12p 
reduces to Shei's single complex kink condensate (14. 5p up to an unimportant constant phase factor. 
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3. Reduction to single plane wave condensate 

In the opposite limit, as 0, the period remains finite, because K 7r/2, and A 2tanh(6'/4). As ^ 0, the 
Weierstrass functions simplify: 

a{x) sin(a;)e^'/6 



-0+—27-J (4-36) 



1 ^ 1 
3 sin^ (x) 

Thus, the complex crystalline condensate ()4.12|) reduces to a single plane wave (chiral spiral) form: 

A ^ msech2(6i/4)e~2»m ta„h^(e/4)x ^4 37^ 

V. SOLUTIONS TO THE BOGOLIUBOV/DE GENNES EQUATION 

A. Spinor solutions 

Remarkably, not only is it possible to solve the NLSE (|2.13p exactly, we can also solve exactly the associated BdG 
equation [i^. In the previous sections we described the spectrum; here we present the explicit spinor solutions and 
express the spectral information in a more compact and useful form. We write the BdG equation (jl.Sp as 

The solutions for the real condensates in Section [1111 are well known For the complex plane wave, A{x) = Ae"'^, 

the solutions are simply chiral rotations of free spinors, and are discussed for example in [44| . For Shei's complex kink 
condensate (|4.5p . the spinor solutions are given in [1, H^. For the complex kink crystal (|4.12p . the two independent 
spinor solutions can be written as 

, / cr(Ax+iK'+ia-ie/4) jj^l-icge /2)+ins(ie /2)]+iem/4+iTr /i 

V 2|dfc/da| I a(Ax+i-K.' +ia+ie/A) i4i£[-iC(i6>/2l+znsfze/2H-i6>m/4-i7r/4 



f^^ d ) = ^V' (5.1) 



(y(Ax+i¥L') (y{ia+ie/i) ^ 



V'4 



A 



2\dk/da\ 



I <y{Ax+iK'-ia-ie/4:) i4^\-ic(ie /2)+ins(ie /2)]+iem/4+i-!v /4 
' a{Ax+iVi')a{ia+ie/i)^ 

, a(Ax+iK.' -ia+ie/4) - i^\-iC{ie /2^+ins{ie l2^]-ierl:J i^-i-n / 4 
\ a(Ax+iK')cr{ia-ie/4:)^ 



e 



(5.2) 

' ^ [iC(ia+ie/4)+zC(ia-ie/4)] 



Here, a is a spectral parameter that characterizes the energy and momentum of these solutions, and k is the momen- 
tum, defined below in (|5.7p . Notice that il'± differ from one another by the sign of a. 

To relate the energy eigenvalue E to the spectral parameter a, we substitute these forms into the BdG equation 
(|5.ip . and make use of the Weierstrass function identity (|11.13p . This identity immediately shows that these are 
indeed solutions, and determines the energy to be 

A 

E{a) = —[iC{ia~ie/A)-iC,{ia + ie/A) + iC{iei2)+ixiii{iei2)] 



V{i6lA)~V{ia) 

where A is given by ()4.13p . The spectral parameter a is defined on the vertical edges, < ia < «K'(i/), K(i/) <ia< 
K(!^) + iK.'{v), of the fundamental rectangle shown in Figure [TUl Only on these edges of the fundamental rectangle 
are the spinor wavefunctions in (|5.2p bounded, and is the energy E{o) real. The vertices of the rectangle correspond 
to the band edges [compare with equation (|4.19p ]: 



Ex = E{ia = 0) = -m 
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FIG. 10: The fundamental rectangle for the spectral parameter ia appearing in the spinor solutions (|5.2p to the BdG equation. 
The bound band is characterized by K(i') < ia < K(!/) + iK'(!/), while the positive and negative energy continua are 
characterized by < ia < iK'(i^). The point a = 6/4 represents E — ±oo, depending on the side of approach. 



Energy Energy 




FIG. 11: Energy E{a) from (|5.3|) as a function of the spectral parameter a, for two different values of 0. The solid [red] curves 
show the continuum energies, while the dashed [blue] curves show the energy in the bound band. The horizontal dashed lines 
denote the band edge energies Ei, E2, E3 and E4, from (15. 4p . and the vertical [red] line gives the asymptote to iJ = =poo at 
a — 9/4, as discussed in the text. These curves are for elliptic parameter v — 0.05. The first plot is for 9 — 2K' = 5.82, so that 
the band is symmetric about the origin. The second plot is for 9 = K.' = 2.91, so the band is asymmetrically offset from the 
origin. 



E2 = E{ia = K) ^m{-l + 2nc^{ie/'i)) 

E3 = E{ia = K + iK') =m{-l + 2n(f{ie/4:)) 

Ei = E{ia = iK') = +m (5.4) 

The right-hand boundary of the fundamental rectangle corresponds to the bound band, while the left-hand boundary 
corresponds to the positive and negative energy continua. The point a — 6/4 is associated with the point at infinity; 
the bottom of the Dirac sea is approached as a ~* 0/4 from below, and the top of the positive energy continuum as 
a 9/ A from above. This is depicted in Figure [TTJ 

To identify the momentum associated with these solutions, we recall that the quasiperiodic winding (|4.14p of the 
condensate in (|4.12p implies that the BdG hamiltonian (|1.7[) is invariant under a period shift, up to a global chiral 
rotation through the winding angle ip: 



H{x + L) = e'^''^H{x) e-*"^^^ 



(5.5) 
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FIG. 12: Momentum as a function of the spectral parameter a, in the positive and negative continua [first plot] and in the 
bound band [second plot]. Note that dk/da is negative in the continuum and positive in the bound band. These plots are 
for v — 0.05 and = K', as in the second figure in Figure [11] The vertical [red] line in the first plot gives the asymptote to 
k = ^oo at a = O/A, as discussed in the text. 



Using the quasiperiodicity properties of the Weierstrass functions (|11.2p . we see that under a period shift, the spinors 
in (|5.2p acquire a chiral rotation and a Bloch phase: 

^±{x + L) = e±^'=^e"^^^V±(a;) (5.6) 

where if is the winding angle defined in (j4.15p . The Bloch momentum k is expressed in terms of the spectral parameter 
a as: 

k{a) = --^ [« C(«a + iS/A) + i ({ia - iO/A) + 2?7a/K] (5.7) 

This is the relativistic version of Bloch's theorem. The momentum is real for a taking values on the vertical edges of 
the fundamental rectangle, and is plotted in Figure [121 Since fc(a) is odd in a, we see that V'lt in (|5.2p are positive 
and negative momentum solutions. 



B. Density of states 

In this Section we present an efficient characterization of the density of states, and show how this relates to the 
resolvent in (|2.12p . In the previous section both the energy E and the momentum k were expressed in terms of 
the spectral parameter a. Now consider the derivatives of the energy and momentum with respect to the spectral 
parameter a. From (|5.3p and (|5.7p we see that 

^ = --[V{ia + ie/A)-V{ia-ie/A)] 
da 2 

dk A 

— = --[V{ia + ie/A)+V{ia-iei'i) + 2i-i/K\ (5.8) 
da 2 

The significance of these expressions becomes clearer once we write the resolvent functions a{E), b{E) and M{E) in 
(|4.20p as functions of a (instead of E) : 

a{a) = —[2V{i0/2)-V{ia + i0/4)-V{ia-i9/A)] (5.9) 
6(a) = A[iC{ia-ie/4)-iC{ia + ie/4:) + iCii0/2)~ins{ie/2)] (5.10) 

■^^"^ " A^V{ia + iO/i) - V{ia - iO/A)] ^^'"^^^ 

Here the upper sign in M is for the continuum states, and the lower sign is for the band. Thus, we recognize dE/da 
as 

dE 



T^ViE^ - m^){E - E2){E - E3) . (5.12) 
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And dk / da can be written as 



dk 

— = -A[P{i6/2)-a{a)/A^ +ri/K\ , (5.13) 
da 



which is negative in the continuum and positive in the band [see Figure [T^]- Consequently, the density of states can 
be expressed as 

dk dk I dE 
dE da da 

2^{E^^m^){E-E2){E-E^) 
When 9 = 2K', which is the GN2 hmit [i.e., with a real condensate], this density of states reduces to 

dE ^(i;2-m2)(£;2_^2(i_j>)) ' (i + V^)^ 

where E and K are the complete elliptic integrals [13, . This density of states (j5.15p is precisely that found for the 
single-band finite-gap Schrodinger system of the GN2 model f?\ . The result in (|5.14p generalizes this density of states 
to the general complex crystalline condensate (|4.12p . 

The density of states (|5.14p has been derived from the energy and momentum of the spinor solutions in (|5.2p . For 
consistency, we compare this with the trace of the resolvent (|2.12p over one period: 



1 r A r+KM 

- / tTDR{x;E)dx = — / 2MiE){aiE) + \A{x)\'')dx 
^ Jl J-k/a 



K 



K/A 

+-K/A 

dx [a{E) + A^ {V{Ax + iK') - V{ie/2))\ 



AN{E) r+-^'^ 



K/A 



[a{E)-A^V{^9/2)-A^r^/K] 



where we have used the integral 



2^{m^~E^){E-E2){E-E3) ' 

+^/^ 2fi 

V{Ax + iK')dx = — f . (5.17) 

This illustrates the consistency of our resolvent ansatz (|2.12p with the spectral properties of the associated BdG 
equation. 



VI. HARTREE-FOCK APPROACH AND THE GAP EQUATION 

So far we have shown that the gap equation (|1.4p can be reduced to the NLSE (j2.13p , and we have found the general 
solution to the NLSE. We now show the self-consistency of this approach by verifying that the gap equation in the 
form p.9p is indeed satisfied. Inserting our ansatz form (|2.12p for the resolvent, the gap equation in this form reads 

A{x) = -2iNghi-E [JV{E) {b{E)A{x) - iA'{x))] (6.1) 

To verify (j6.ip . we require information about the single-particle fermionic spectrum of the BdG Hamiltonian for the 
condensate A{x) that we have found. This provides a complementary approach to finding a self-consistent condensate, 
by solving the associated relativistic Hartree-Fock problem. 

A. Solving the gap equation from the resolvent 

We first need to show that the coefficient of A'(a;) on the right-hand-side of (|6.ip vanishes once the energy trace is 
taken: 

tiE WiE)] = (6.2) 
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FIG. 13: The occupation of single-particle fermionic states for two different values of 6. In general, to satisfy the gap equation 
the condition (|6.5I) relates the condensate parameter 6 to the filling fraction of the bound band. In the special case 6 — 2K', 
this filling fraction is 1/2, and this corresponds o the real condensate of the discrete-chiral GN2 model, for which the spectrum 
is symmetric about 0. Thus, in this case, all the negative energy states are filled. 



This computation is greatly simplified by converting it into an integral over the spectral parameter a. Recalling (|5.12p 
and the normalization factor in (|4.20p we have 

^E) = - (6.3) 

We assume that the negative energy continuum is fully occupied, and that the band inside the gap is partially occupied, 
by a fraction ^. Then we can express (16.21) as 

~ 7 If da da> 

= ^{e/i^W (6.4) 

Thus, 9 must satisfy the condition 

— — = ^ = fiUing fraction (6.5) 



To appreciate this condition, we consider some special limiting cases. In the real [GN2] limit, where 9 = 2K', this 
corresponds to filling the band half-way. Since the band is symmetrically placed about zero energy, this corresponds 
precisely to filling all the negative energy states, as in the conventional GN2 model analysis 0, Q. For the single kink, 
which is obtained in the infinite period limit, the bound state is half occupied (28l. Is^. 

For general 9, in the infinite period limit, where K' —^ 7r/2, we recover Shei's condition ([4lT|l that 9/27r is the 
fractional filling of the bound state in the gap. Note that in this case, the fractional filling of this state is conventionally 
interpreted [1, jSj] in terms of the filling of the level by a fixed fraction of the fermion flavors, with this fraction 
kept fixed in the infinite N limit. At finite period, it is thermodynamically more natural to consider a fraction of the 
band being occupied by all flavors, as above. As the band contracts to a single bound state in the infinite period 
limit, we smoothly tend to the situation of having the level partially filled, while maintaining consistency with the 
gap equation for any period. 

The other condition required for the gap equation (|6.ip to be satisfied is that 

- 2iNghvE [M{E) b{E)] = 1 . (6.6) 
This is the standard vacuum gap equation for the renormalization of the coupling 0, H, 0) IE] ■ 



B. Solving the gap equation from the spinor solutions 



Given the spinor solutions of the BdG equation, we can reconstruct the expectation values (V'(a;)V'(a;)) and 
— {■ip{x)i'^^tp{x)), to verify that they correspond to the real and imaginary parts of condensate A{x). From (|5.2p . 
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we find 



.-2^rJ,a^Ax[-^C{^e/2)+^nsi^e/2)] '^i^^ + + " ^g/4)g(Aa: - iK' ~ ia - i9 / 4) 



\dk/da\ a{Ax + iK')a{-ia - i0/A)a{ia - i0/A)a{Ax - iK') 

= Tirm [C{Ax + zK') + C{ia - ie/4) - CiAx + iK' - 19/2) - C(m + iO/A)] A(.t) 
\dk/aa\ 

where the subscript k on tpk emphasizes that these are spinor solutions of a given momentum k. In deriving (j6.7p we 
have used the product identity (|11.13p for the Weierstrass sigma fimction. 

Thus, for a given momentum fc, ipkix)ipk{x) — i ipkix)ij^ipkix) has the same x dependence as the upper off-diagonal 
entry of the resolvent, which appears in the gap equation (|6.ip . To compute the expectation value we trace over the 
N flavors and integrate over momentum: 

/dk - 
— [i>k{x)%l;k{x) - iil)k{x)i^^ilJk{x)\ 

NT, f dk 



2ttA 



j dasgn (^^^ [b{a)A{x) ~ zA'(x)] (6.8) 



Thus we have precisely the same integrals as were evaluated in the previous section to show that the gap equation is 
satisfied. The vanishing of the coefficient of A'(x) leads to the same relation (16. 5p between the filling fraction and 6*, 
while the the coefficient of A(x) is again the vacuum gap equation (16. 6|) expressing the renormalization of the coupling 
constant g. 

It is interesting to note that while the condensate A (a;) is spatially inhomogeneous, the expectation value of the 
charge density is spatially uniform. From the explicit spinor solutions in (j5.2p we find: 

For a given single particle state of momentum fc, the charge density '\\)\{x)'il)k{x) is spatially inhomogeneous, but this 
inhomogeneity is washed out by integrating over fc, so that the expectation value is uniform: 

-^(V'^(x)V'(a;)) = (6.10) 
dx 

Indeed, for a given fc, the charge density '\\)\,(x)'il)k(x) has spatial inhomogeneity determined by |A(x)p, the magnitude 
squared of the condensate. But the coefficient of this x-dependent term in (16. 9|) is precisely the same as the coefficient 
of the A'(a;) term in (|6.7p . We just saw above that this term vanishes [when the energy trace is taken], in order to 
satisfy the gap equation. Thus, the very same condition (j6.5p that ensures that the gap equation is satisfied also 
shows that l^^'^ {x)il){x)) is uniform. 

There is another, more physical, way to understand this fact {4^ . It expresses the conservation of axial charge in 
the NJL2 model. Recall that the NJL2 system has two conserved currents: the first is the charge current — tj^j'^ip, 
and the second is the axial current = ip^f^^^ip. In 1 + 1 dimensions these are related by = e^'^ j^- For static 
but spatially inhomogeneous condensates, the charge current conservation is automatically satisfied. Axial current 
conservation is more interesting. For a static condensate, 

so axial charge conservation requires a uniform charge expectation value, as was found above (|6.10p from the single- 
particle spectral properties. Another way to see this is to observe that we can express the expectation value of the 
axial current in terms of the resolvent 

(j-^) = tTE.D{fri''R{x;E)) (6.12) 
Taking the spatial derivative and using the Eilenberger equation (|2.6p . we find 

^{jD = trE^D{dMx-E)) 

= 2G{x){i){x)i'fi}(x))-'l'K{x){ii{x)^{x)) (6.13) 



(^)> = ^M'\x)^^(x)) (6.11) 
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Thus, once again, the gap equation ensures that axial current conservation is satisfied. Interestingly, it is encoded 
into the Eilenberger equation (|2.6p . 

On the other hand, in the GN2 model, which has just a discrete chiral symmetry, there is no axial current conser- 
vation and tt{x) = 0. Then the Eilenberger equation (|2.6p again expresses the correct relation 

d^{j^) = 2aix)mx)ij^ijix)) (6.14) 

and {ip{x)'ip{x)) , (j{x) and (7/>(a;)i7^^(a;))are all inhomogeneous. 



VII. ALL-ORDERS GINZBURG-LANDAU EXPANSION 



In this Section we present another perspective on our solution, in order to illustrate what aspects might possibly be 
extended to higher dimensional inhomogeneities. While a direct application of the resolvent approach, or the inverse 
scattering approach, to higher dimensions is problematic, one approach that can be straightforwardly generalized to 
higher dimensions consists of the Ginzburg-Landau expansion of the effective action. In this approach, one expands the 
free energy corresponding to the effective action (|1.3p in powers of the static condensate A [i.e, the order parameter] 
and its spatial derivatives. Using the relation (|2.2p between the spectral function p{E) and the trace of the resolvent, 
we write the free energy density (per flavor) as 



£h-coix) = ^ / dElm [tiD R{x; E + ie)] In (l 



(7.1) 



Thus, knowing the resolvent means we know the free energy. Usually, it is impossible to evaluate this local density 
of states exactly. However, an asymptotic expansion of R{x] E) can always be obtained from a Laplace transform of 
the heat kernel expansion, leading to 



i?(x;i?) = -)^-^ 

n=0 



(7.2) 



We stress that such an expansion can be derived in any dimension, not just in 1 + 1 dimensions. This is because 
the heat kernel expansion is known in any dimension. However, in 1 + 1 dimensions we now have the luxury of 
having found an exact solution. In this Section we compare this exact solution with the Ginzburg-Landau expansion, 
and show in explicit detail how the gap equation for the inhomogeneous condensate is solved order-by-order in this 
Ginzburg-Landau expansion. The way this works is surprisingly simple and elegant, as shown below. 

In 1-f 1 dimensions, including the Hartree-Fock double-counting correction, and renormalizing, leads to the standard 
Ginzburg-Landau expansion of the effective Lagrangian, the low order terms of which are 



C 



GL - ao + «2| AP + «3lm [A(A')1 + a4 [|A|4 + |Af ] + a^lm [(A" - 3|ApA) (A')*] 
+a6 [2|A|6 + SlAplAf + 2Re ((A')'(A*)2) + |A'f ] + . . . 



(7.3) 



The low order terms are relatively simple, but high orders rapidly become cumbersome. The coefficients Q!„(T, /x) are 
known functions of temperature and chemical potential. For example, in 1 -t- 1 dimensions (46l . |47| 



ao 
as 



2tt 



6 



1 

2^ 



ln(47rT) + Re^ [l+i^ 
2 27r 



2^TT^T ^ I 2 27r 



1 



267^3^2 



ReV/ 



2+^2. 



1 



2i27r53r4 



Re 



'2^ 



(7.4) 
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Here tp^'^'^ denotes the k*^ derivative of the Euler digamma function V'(-z) = d\nr{z)/dz. 

In higher dimensions, an analogous Ginzburg-Landau expansion can be derived. The form of the q;„(T, /i) is 
different, but known, and the form of the spatial terms r„(a;) is different, but computable. In this section we consider 
how the NJL2 model gap equation for inhomogeneous condensates looks in terms of the Ginzburg-Landau expansion. 
We first present a simple recursive way to generate the expansion to all orders, and then we show how the gap equation 
is solved order by order. A closely related expansion, the derivative [or gradient] expansion, which is an expansion 
just in powers of derivatives, but including all orders in powers of the condensate, is considered in |36i] . 



A. Ginzburg-Landau expansion in 1 + 1 dimension 

The Ginzburg-Landau expansion of the effective action follows, by (|2.8p and (|2.2p . from an asymptotic expansion 
(|7.2p of the resolvent R{x; E). The Eilenberger equation (|2.6p generates such an expansion recursively. It is more 
convenient to define (as in [36|) 

gix;E) = R{x;E)a3 (7.5) 
In terms of g, the Eilenberger equation (12. 6p reads 



g' ^iE[a3,g]-i[J,g] , J = ( ^] (7.6) 



We now define an asymptotic expansion for g 



00 



71 — ^ ' 

The factor ~i comes from the large E expansion of J\f{E), as is already clear from the constant A case in p. lip , 
which also tells us that go — <^3- The Eilenberger equation implies the simple recursion formula for the gn{x): 

g'„{x) = i[a3,g,i+i{x)] - i[J (x) , gn{x)] (7.8) 

The determinant condition det g — \ fixes the c„ in terms of the c?„ : 

Co = 1 

ci = 

C3 = ^ {d2d*i + did*2) 

C4 = ^\\d,\^ + \\d2\^ + \{d3dl+d,dl) 

C5 = -\\di\^ {d2dl+did*2) + ]^{d2dl + d3dl2) + ]^{didl+didl) 

'■■ (7.9) 

So, given ga — 03,, we learn from (|7.9p that ci = 0. Then the Eilenberger recursion equation (j7.8p determines c?i = A, 
so that gi{x) — J{x). Next, knowing di, we learn from (17. 9p that C2 — ||Ap, and from the Eilenberger recursion 
equation (17. 8p . we find that c?2 = — Iterating this procedure we find: 

gi{x) = J{x) 

1 /lAP iA' 



92{x) 
g4{x) 



2 \iA'* -|Ai2 

_1 /i(A'A* - A'*A) -(A"-2|A|2A) 
4 lv(A"*-2|A|2A*) -i(A'A* - A'*A) 

1 /3|A|4 + 3|Af - (|A|2)" -t (A'" - 6|A|2A') 
S\ -z(A*"'-6|A|2A*') -3|A|4-3|Af + (|A|2)' 



(7.10) 
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The spectral function is expressed in terms of the trace of the resolvent: 

trR{x■E)^-^Y^'-^ (7.11) 

n=0 

Then the Ginzburg-Landau expansion (|7.3p is obtained by inserting (|7.1ip into (|7.ip and performing the energy 
integrals. Renormalization affects the first two terms, and the others follow from the integrals {n > 2) |4^]: 



2)/2^„-2 .^^ ( i + i Hii ) , n even 



^/ d£;im ^ ) log ( 1 + e-'5(^-^n ^ <^ ' ' ' ' ^ ^ 



B. Gap equation to all orders in Ginzburg-Landau expansion 

We now show how the gap equation (|1.4p is solved order-by-order, to all orders in the Ginzburg-Landau expansion 
(I7.12p . There are two ways to see this. First, even though the c„(a;) at high order n become extremely complicated 
combinations of powers and derivatives of A and A*, when we evaluate Cn{x) on a solution of the NLSE (j2.13p . they 
dramatically simplify and, for each n, reduce to something linear in |A(a;)p: 

Kix)]NLSE = l3n\Mx)\^ +ln , (7.13) 

for some constants /3„ and 7„. For example, for 02(2;) = |A(a;)p it is obvious, and for C3(x) it follows immediately 
from (gSlI): 

C3(x) = -1(A'A*- A'*A) 

= -i(6- 2^;) I Ap - constant (7.14) 

For C4, we note that the NLSE ([2^3]) implies that (|A|2)" = 4|A|4 + 2|A'|2 -t- 4 (a - P/4 ~ E'^) |A|2 + constant, so 
that 



C4(x) ^ ^(|A|4 + |A'p-i(|AH") 



= i (-|A|4 + |A'|2 - 4 (a - 6V4 - |A|2 + constant) (7.15) 
8 

This expression is linear in |Ap because we recall from (I^Ti)) and (ji?^ that (-|A|'' + |A'p) is linear in jAp; this 
ensures that AeiR{x\E) = —j is satisfied, and shows that 04(2;) reduces to the form in (|7.13p . Indeed, (|7.13p holds 
to all orders in the Ginzburg-Landau expansion - the diagonal component of the Eilenberger equation (|2.6p is the 
generator of this infinite sequence of relations. 

Another related fact is that when we seek stationary points of the Ginzburg-Landau effective action with respect 
to variation of the condensate, we generate a sequence of equations for the condensate. At higher order n, these 
equations become more and more complicated, involving higher derivatives and nonlinearities. But, when evaluated 
on a solution to the NLSE, these variations take a much simpler form, for all n: 

dx Cn{x) ^ 5n /^{x) + tn ^' [x) (7.16) 



5A*(x) 



for some constants 5n and tn- It is straightforward but instructive to verify this remarkable property for low orders. 
To all orders, this infinite sequence of identities is generated by the off-diagonal entries of the Eilenberger equation 
(|2.6p . so the off-diagonal entry of R{x\ E) in (|2.12p explains why (|7.16p holds at any order. To solve the gap equation, 
the parameters of the solution [i.e., the parameter 6] must be adjusted so that the coefhcient of A'(x) vanishes, as 
explained in Section IVTl 

These observations may be a useful guide to seeking solutions to the gap equation in higher dimensions where 
neither inverse scattering nor the resolvent approach is directly applicable for periodic condensates, but where the 
Ginzburg-Landau expansion approach is available. 
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VIII. CONCLUSIONS 



To conclude, we have solved the inhomogeneous gap equation (|1.4p in the NJL2 model, to obtain a self-consistent 
condensate (I4.12p that is crystalline in form. It is inherently complex, unlike the condensate in the discrete-chiral GN2 
model where the condensate is real. The complex crystal has a periodic amplitude, while its phase rotates though a 
certain angle over each period. This complex crystalline condensate contains all other previously known solutions to 
the gap equation as special cases. In addition to finding the exact condensate, we have presented the exact solution to 
the associated Bogoliubov-de Gennes equation (jl.Sp which gives the spectrum of fermions quantized in the presence 
of such a condensate. The spinor wavefunctions and single-particle energy spectrum are found in closed form, and we 
confirmed the consistency of the results by solving the gap equation also by Hartree-Fock. The key technical idea in 
our approach is to use the form of the gap equation to motivate an ansatz (I2.12p for the Gorkov resolvent. For any 
condensate, the resolvent must satisfy the Eilenberger-Dikii equation (12. 6p , and so this reduces the problem to solving 
the nonlinear Schrodinger equation (|2.13p . Given this solution, and in particular the exact density of states (|5.14p . 
we are now in a position to study exactly the phase diagram of the NJL2 model, extending the Ginzburg-Landau 
analysis in For example, once the free energy is minimized on our exact gap equation solution, the filling factor 
is then determined as a function of T and /i. This will be addressed in future work. 

An important possible extension of this work would be to the massive NJL2 model, in which a bare fermion mass 
is included in the original Lagrangian (jl.ip . This mass term explicitly breaks the continuous chiral symmetry, and 
approximate methods have found a rich structure in the associated phase diagram [i^, [4^] . The resolvent approach 
presented here is easily extended to the massive case [yielding the real kink-antkink crystal condensate of the massive 
GN2 system] only when the condensate is real. This is deeply related to the integrability properties of the classical 
equations of motion of these systems (soj . More work is needed to understand analytically the situation relevant 
for the massive NJL2 system, where the condensate is complex. To date, no exact self-consistent inhomogeneous 
condensate solving the gap equation has been found for the massive NJL2 model, except for simple embeddings of 
the real solutions 51]. Another possible extension is to study other 1-1-1 dimensional models, such as the Schwinger 
model [52, |531. 

The most interesting extension would be to try to extend some of these ideas to higher dimensions, for example 
for higher dimensional Gross- Neveu or Nambu- Jona Lasinio models 5J] , or more ambitiously to search for crystalline 
condensates in QCD or QCD models [H, [H, [13, [H, [s^ . Here, the resolvent approach is not directly applicable, but 
the lessons from the way in which the Ginzburg-Landau expansion works (discussed in Section IVIip may be useful. 
A numerical approach currently being developed (60j . that does not rely on special features of 1 -|- 1 dimensions, 
is a numerical evaluation of the free energy, using the worldline Monte Carlo approach of Gies and Langfeld [6]| . 
Other new lattice methods have been developed recently |62, |63|] to study four-fermion interactions, and it would be 
interesting to extend these to higher dimensions. 
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X. APPENDIX A: RESOLVENT OF A DIRAC SYSTEM, AND THE EILENBERGER EQUATION 

In this appendix we summarize the derivation of the Dikii-Eilenberger equation (j2.6p for the diagonal resolvent of 
a Dirac system. The interested reader is urged to consult also references [SJ, |35j , and Appendix B of the first paper 
in [3^. Here we sketch the key features of the argument, for the sake of being more self-contained. 

The important idea is simple. Recall that in one dimension, for a Schrodinger-like Sturm-Liouville operator, it is well 
known that the Green's function can be written in terms of a product of two independent solutions, normalized by their 
Wronskian [38| . An analogous construction exists for a one dimensional Dirac operator 35, 36J. Furthermore, in 
the Sturm-Liouville case, the coincident-point limit R{x, x; E) satisfies a differential equation, known as the Gel'fand- 
Dik'ii equation (for an excellent review, see Q), just by virtue of being written as a product of solutions to the 
original differential equation [3^. Likewise, for a Dirac system, the coincident-point limit R{x,x;E) also satisfies an 
equation, just by virtue of being expressed in terms of solutions to the original differential equation. This equation 
is the Dik'ii-Eilenberger equation (|2.6p . The technical difference from the Gel'fand-Dik'ii equation arises because the 
Dirac operator is first- order in derivatives and because it is a 2 x 2 matrix operator. 

So, consider two independent spinor solutions ^pl,2 of the BdG equation Hijj = Eip in (|1.8p . where H is the 2x2 
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matrix first-order differential operator in (|1.7p . Then we can write the resolvent (also a 2 x 2 matrix) as 

The row- vector functions R{y) and L{y) are determined by demanding that R{x,y;E) satisfy 

{H ~ E)R{x,y;E)^5{x-y) (10.2) 
Integrating this first-order differential equation across the point x = y, we learn that 

where the Wronskian W is the scalar function W{ipi,ip2) = i4'i'^2'>p2- The coincident-point limit x = y is given by 
an averaged limit as [see p.7p ]: 

Rix; E)^^ {Mx)^p^ix) + Mx)i^Iix)) cTi (10.4) 



Now, express the BdG equation (jl.Sp as 

^' = ^(f. -_e)^ ' ^'^ = #^(_1 (10.5) 

Then, differentiating (|10.4p once with respect to x, using the fact that each of tpi and ip2 satisfies the BdG equation 
p.8p in the form p0.5p . we arrive immediately at the Dik'ii-Eilenberger equation (|2.6p . Furthermore, if we write 
-01 = (mi,i;i)-^, and 7/12 = (^2,^2)'^, then p0.4p says that 

R{x;E)^- ^ (u,V2 + U2V^ 2u,U2 \ 

2i{uiV2 - U2V1) \ 2viV2 U1V2+U2V1J 

Thus, the algebraic conditions (|2.4p and (|2.5p follow. The hermiticity condition (|2.3p is less immediately obvious, 
because this depends on E and W. However, it is clear from the definition of the resolvent that R{x; E) = (x\l/(H — 
E)\x) is hermitean for real E. The appropriate hermiticity properties are studied in more detail in [341 l35l l3q. and 
can also be seen from the result for the constant A case (|2.1ip . and the ie condition used to define the spectral 
function ([2?2l) . 

XI. APPENDIX B: SOME USEFUL PROPERTIES OF ELLIPTIC FUNCTIONS 

In this appendix we collect some basic facts and nontrivial identities for Weierstrass elliptic functions that are 
used repeatedly in this paper, in order to make the paper more self-contained. These functions play a special role 
because the self-consistent condensate A (a;) and the spinor solutions ip{x) to the Bogoliubov-de Gennes equation 
are all expressed in terms of these functions. There are many good books on elliptic functions. Excellent classical 
references are [13, HI] . We also found [s^ H^] to be particularly useful. Very roughly speaking, the Weierstrass elliptic 
functions are doubly-periodic extensions of standard trigonometric functions: 

sin(z) ^ a{z) 

cot(z) = — Insin(z) (^{z) — —\na(z) 
dz dz 

^ - A cot(z) ^ Viz)^ --few (11.1) 
sm (z) dz dz 

The trigonometric functions have periodicity properties along the real z axis, but the Weierstrass functions are doubly 
(quasi-)periodic. They are specified by the real and imaginary (half-)periods, lui and W3. It is standard to define also 
UJ2 by + C1J2 + ^3 = 0. Then the Weierstrass sigma function is quasi-periodic under shifts by 2^^: 

a(z + 2c^,) = -e"'*^"+"'V(z) , z = 1,2,3 . (11.2) 
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Here rji = C(wi). In general, uJi and W3 define a fundamental parallelogram characterizing the doubly-periodic nature 
of the Weierstrass functions. We choose a fundamental rectangle, with cui — K(zy) and u;^ ~ iK.' = i K(l — i/). The 
periods are then parametrized by the elliptic parameter 1/ that takes values in [0, 1]. At the limits v—O and 1^ = 1, 
the Jacobi elliptic function reduce to trigonometric and hyperbolic functions respectively. Physically, they interpolate 
between kink-like solutions and sinusoidal ones. 

The Weierstrass functions arc then related to the Jacobi elliptic functions as follows. We define 

..W = e--^(^ (11.3) 

and Jacobi's elliptic functions can be constructed from their ratios, 

sn = —f-T ; cn(z) = — -- ; dn(z) = — -- (11.4) 
(y^yz) a-i{z) a-i{z) 

where the Jacobi functions have elliptic parameter v and the Weierstrass functions have periods lo\ — K(z/) and 
Weierstrass's zeta function is defined as the logarithmic derivative of the sigma function: 

C(.) = ^lnK.))^^ (11.5) 

From (|11.2p it is clear that C is also quasi periodic 

C(z + 2cj,) - 2r;, + C(z) ; i = l,2,3 . (11.6) 
Here r\i is given by the zeta function evaluated on the periods z = Wi: 

aio,)=r,, (11.7) 

Finally, the Weierstrass V function is defined as: 

n.,--^ ,u.a, 

V is doubly periodic, with periods 2wi, 2^3: 

Viz + 2uji) ^Viz) ; i = l,2,3 . (11.9) 

Another important property of V is that it satisfies the differential equation 

V'^iz) = AV^{z)-g2P{z)-g3 

= 4{V{z)~e,){V{z)-e2){V{z)-es) (11.10) 

This equation is the one (|4.22p satisfied by the amplitude squared of the condensate, following from the nonlinear 
Schrodinger equation. The constants 52 and 173 in (|11.10[) are known as the invariants and they are parameters 
depending on the periods. Similarly for the e^, with £1 + 62 + ^3 = 0, and 52 = — 4(eie2 + 6263 -I- 6163) and 
gs = 4616263. With our choice of periods, these can be related to the Jacobi elliptic parameter v as: 

6i = i(2-i.) ; 62 = i(2z/~l) ; ei = -^{1 + 1^) . (11.11) 

Just as the trigonometric and hyperbolic functions satisfy addition and product formulae, so too do the elliptic 
functions. Indeed, these elliptic identities generate all others as special cases. Here we list some of the important 
addition formulas for a, C and V that we have used throughout the paper. References [3^ [4^ are particularly good 
concerning these identities. 

^^^^i^fr^— ^(-)+w (11.12) 

(7^(u)0"^(w) 



(11.13) 



[C{u + v)~ C(") - Civ)] = Viu + v)+ V{u) + V{v) 
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(11.14) 



V'iv) 



T{v)-V{u) 



(11.15) 
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